RESCEU-34/11 



Evolution of FLRW spacetime after the birth of a cosmic string 

Matthew Lake^ and Teruaki Suyama^ 

Research Center for the Early Universe (RESCEU), 
Graduate School of Science, The University of Tokyo, Tokyo 113-0033, Japan 

(Dated: May 2, 2012) 

We consider the evolution of an initially FLRW universe after the formation of a long, straight, 
cosmic string with arbitrary tension and mass per unit length. The birth of the string sources scalar 
and tensor-type perturbations in the background metric and both density and velocity perturbations 
in the background fluid, which compensate for the string mass and maintain energy conservation. 
The former generate the deficit angle within the light cone of the string and a gravitational shock 
front at the cosmological horizon, whereas the latter are confined within the sound cone. We study 
the properties of the metric within each region of the resulting spacetime and give the explicit 
04 . coordinate transformations which demonstrate non-violation of causality. This paper generalizes 

the work of previous studies for the Nambu-Goto string. 
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J>^- I. INTRODUCTION 

^[ 

Cosmic strings are linear concentrations of energy which may have formed during symmetry-breaking phase tran- 
sitions in the early Universe Though the string width is determined by the inverse of the symmetry-breaking 
energy scale, this is smah compared to cosmological distances, and they may be approximated as one-dimensional 
objects for many purposes @, Q- Specifically, strings may have been produced at the epoch of electro- weak symmctry- 
^ ■ breaking , or the GUT scale, but their formation is also a generic feature of the phase transitions inherent in many 
extensions of the standard model, (c.f. @ and references therein). In field theory they are a type of topological 
5-H . defect, analogous to the magnetic flux tubes and other vortex-type defects created in certain condensed matter sys- 
tems |10l - [l§ |. and are produced via the Kibble mechanism [l^, if the vacuum manifold {M.) possesses a nontrivial 
first homotopy group (e.g. ■ni{M) = Z, in which each element, a nonzero integer, corresponds to an allowed winding 
CNJ ' number for the string vortex cross-section). In recent years, the formation of "cosmic", i.e. horizon-sized, fundamen- 
^ , tal strings (i^-strings) and one-dimensional ZJ-branes (D-strings) has also been extensively studied in string theory 

■ [20l - l30| . particularly in the context of brane inflation (see fsTI - fSq for reviews), in which such defects can be copiously 
I produced (s^-fl^- However, regardless of the precise details of individual models, the main phenomenological, and 

^ij^ . observationally relevant, parameters that characterize the string are the energy of the per unit length, which we denote 
U, and the tension, n = qU . ^ The Nambu-Goto (NG) action represents the simplest model and enforcing Lorentz 
. symmetry along the string corresponds to setting g = 1. If the NG string is wiggly at microscopic scales, it may be 
approximated macroscopically by a simpler embedding with an effective mass per unit length which differs from the 
intrinsic tension. Thus in the coarse-grained limit, small-scale structure, even on pure NG strings, can lead to a value 
of q different from un ity [49l [50|. For current carrying strings, the mass-density and tension may differ, as originally 
_ ^ _ shown by Witten [5l|, [53| , so that g 7^ 1 in general. 

■ The gravitational field around a static, straight string exhibits an interesting feature not seen in fields surrounding 
\^ ' spherically-symmetric distributions of matter. In a perturbative analysis of a string in the wire approximation (to 
5^ . first order in G/x), Vilenkin [H^l showed that test particles within the string light-cone experience a gravitational 

acceleration along the radial direction equal to 2GU{q— l)/r, where r is a distance from the string, and that the sur- 
rounding space has a conical structure with deficit angle 47r(l -I- q)GU. For the NG string therefore, the gravitational 
force vanishes and the spacetime becomes locally Minkowski. These results were found to hold in a nonperturbative 
analysis (i.e. to all orders in Gyu) , independently, by Gott [s^ and Hiscock [s^] , even for a string of finite width. One 
observational consequence is the formation of double images from bright sources lying behind the string in relation 
to an observer. The images have exactly the same shape and brightness and for NG strings their angular separation 



This is not strictly true in the case of models containing compact extra-dimensions in which the intercommuting probability for strings 
which "cross" in the noncompact space can be much less than unity (P << 1). In this case the dynamics of string networks and 
their implications for cosmological observations can be substantially altered [sol. |37|440|| . Extra-dimensional scenarios, especially those 
motivated by, or directly embedded into, string theory also allow for the formation of cosmic "necklaces" of various kinds, which 
introduce significant alterations to the network dynamics along with additional phenomenological parameters, such as the bead mass 
|41|446|| . The analogues of a superstring necklaces in field theory are series of monopoles connected by string-segments, which may also 
be formed in some symmetry-breaking models l47l |48|| . The authors hope to address the important question of the inborn metric for a 
cosmic necklace in a future article. 
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on the sky is of the order of the deficit angle , thoug h for g >> 1 it may be substantially smaller [56|, |57[. This feature 
has been used to search for long strings [Hj, [58l - l66l | and loops [67l - [69| via astronomical observations. Bounds on phe- 
nomenological parameters such as the string tension and number density have also been obtained using the predicted 
effects of strings on CMB anisotropy [70l - [89i. the expected emission of high-energy rays from cusps and kinks [90l493| 



and possible signatures in the 21cm line 9J|. Interestingly, it has recently been proposed that cusp-emission from 
superconducting strings may be responsible for anomalous gamma-ray burst observations [qsI - IqtI I and there remains 
ongo ing controversy on that point (c.f. 98l - ll06 | and references therein). With the advent of the Planck experiment 



[l07| , it is also hoped that the polarization of CMB S-modes caused by s tring netw orks in the early universe could 
confirm, or rule out, the existence of GUT -scale str ings in the near future 108l4ll2j |. Finally, gravitational radiation 



from strings ha s be en extensi vely studied jll3l - lll7 and future gravitational wave detections from experiments such 



as LISA/NGO |118l | or LIGO |119l | may place even stronger bounds on string parameters. 

It is interesting to consider how the conical structure of the spacetime surrounding the string evolves from the initial 
Friedmann-Lemaitre-Robertson- Walker (FLRW) background, which exi sts t j efore the symmetry-breaking epoch, and 
a detailed analysis of this problem was originally given by Magueijo jl20l Il2l| for the case of an NG string. He 
introduced a phenomenological string energy-momentum tensor which is localized at the origin and appears only after 
the symmetry-breaking time. The generation of string mass-energy is compensated by a decrease in the mass-energy 
of the perfect fluid that dominates the Universe, so that conservation of energy and momentum are not violated. ^ 
Inclusion of the background fluid is therefore compulsory and the problem of string generation must be considered 
on an (initially) FLRW background. In his analysis, Magueijo obtained solutions of the linearized Einstein equations 
on the FLRW background metric by assuming the dimensionless parameter GU to be small and constant. This is 
consistent with known observational bounds on the string tension in various cosmological models, which give an upper 
limit of GU < 0(10"^) for field -the oretic strings [2|. [s^ - ls?! and 10~^^ < GU < 10~^ for the warped tension of cosmic 
F/_D-strings [s^ [s^ (see also |l25| for a review of observational constraints up to 1990). He showed that, at late 
times and in the re gion close to the string, the spacetime takes the form of that surrounding an eternal string having 
deficit angle SttGU |53l - [55| and that a gravitational shock propagates along the lightcone, with the spacetime outside 
the horizon remaining unperturb ed s o that the FLRW metric remains valid, though detailed calculations confirming 



this picture remain unpublished |121 



In this paper, we reconsider this issue and also extend the original analysis by considering the more general case of 
q ^ 1. This allows us to include the physically interesting class of "conical" spacetimes which are not locally fiat and 
to consider the gra vitat ional acceleration of test particles in the vicinity of a noneternal string. We follow the same 
approach taken in |l20l | to evaluate the metric perturbations induced by string formation and find that the general 
picture remains qualitatively similar. However, as shown in detail in Sec. IV, our solutions do not exactly match 
those obtained previously, even for q = 1. We believe that this i s be cause of the choice of boundary condition for the 
Poisson equation governing the gravitational potential given in |l20l | is not consistent with trace part of the Einstein 
equations. We also provide concrete coordinate transformations both for inside and outside the sound horizon of 
the perfect fluid, as well as in the region between the sound horizon and the edge of the light cone, which result 
in an FLRW metric beyond the causal horizon and the metric f or an eternal string deep inside the sound horizon, 
respectively, in accordance with previous results [ssl - IHsl . Il20l Il2lj . 

The structure of the paper is as follows: In Sec. II we define the perturbation variables and decompose them into 
scalar and tensor components. In Sec. HI we derive the master equations for each typ e, wh ile analytic solutions are 
given in Sec. IV together with a comparison of the corresponding solutions obtained in |l2C| . Section V is devoted to 
studying the properties of the resulting class of spacetimes and a brief summary of our results and discussion of the 
prospects for future work is given in the conclusion. Sec. VI. 

II. DEFINITION OF THE PERTURBATION VARIABLES 
A. Background dynamics 

We assume that the background spacetime is an FLRW universe with metric 

ds'^ = a{i]Y {-djf + jijdx'dx') , (1) 



^ In previ ous analy ses this was not always the case, though particle production by cosmic strings was studied in detail in the following 
sources, [121^121 . 
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where 7^ = diag(l, 1, r^) is the spatial part in cylindrical coordinates and rj is the conformal time coordinate, defined 
via a(r])dr] = dt. We also assume that the universe expands due to the presence of a perfect fluid with energy- 
momentum tensor 

™Tt = (p+PK«. + pn- (2) 

where p is the energy density, P is the pressure and, from the normalization condition u'^u^ = — 1, we have that 

1 



a 



0,0,0 , = (-a, 0,0,0). (3) 



The Friedmann equation and continuity equation are given, respectively, by 

n^^^pa^ p + 3{l + cl)np = 0, (4) 
where ' indicates the derivative with respect to conformal time and from these two equations, we find that 

""=0^1 H^-i(l + 3c^)H^ (5) 

B. Metric perturbation 

The scalar part of the perturbed metric may be diagonalized using gauge degrees of freedom and the full metric 
can be written in the form, 



- (1 + 2<P)dif + (1 + 2'<il)"f,jdx'dx^ + hijdx'dx^ 



(6) 



The gauge invariant tensor-type perturbations hij satisfy the transverse and traceless conditions, 

D.h'j = h\ = 0, (7) 

where Di is the covariant derivative defined for the three-dimensional metric 7^ . We do not consider vector-type 
perturbations since they are not generated by the situation we are interested in (see Subsec. Ill Cp . The nonzero 
components of the perturbed Einstein tensor for the scalar-type perturbations are given by 

SG°o = \(3n^^ + D^'f-3n'^y (8) 
2 



<5G", = -2 A - n^) , (9) 



^G', = ^ {d'D, U^^D') + -i^<5Vn„ (10) 



where we define Sg and llg according to 

I]g = -$-^', Hg = 2D^(4> + *) - 6 ( -H'^ - 2- ) $-F6-H ($ - 21') - e^-. (11) 



For the tensor-type perturbations, it transpires that only diagonal components are excited by the string energy- 
momentum tensor, so that we may set K^^ = h^g — h^g — without loss of generality. The tranverse and traceless 
conditions also allow us to write and h^g solely in terms of h^'^, so that h^j, is the only independent component; 

h% = h\ + r—h\, h\ = -2h\ - r—h\. (12) 
or or 

Rewriting h^^ = the relevant component of the perturbed Einstein tensor is the r — r component, given by 

5G\- = A + 2H/^ - F" - -F'] . (13) 
za^ V r I 
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C. Energy-momentum tensor of a cosmic string 

As mentioned previously, in field theories, the string width is roughly equal to the inverse of the symmetry-breaking 
energy scale which, in realistic models, is likely to be above the electroweak scale For cosmic superstrings it is 
expected to be of the order of the fundamental string scale which, although unknown, may be as small as the Planck 
length [2l|. Since both of these length-scales are typically much smaller than the Hubble horizon, the finite width of 
the string does not play a cosmologically significant role as long as we are concerned with purely gravitational effects. 
The finite duration of the string formation as a result of the phase transition is also much shorter than the Hubble 
time and is not significant on cosmological timescales. As considering these effects would complicate the analysis 
of our problem immensely, we avoid such complexity by assuming that the string width is zero (adopting the wire 
approximation), and that string formation occurs instantly. 

With this approximation, the energy-momentum tensor of the cosmic string can be written as 

U6(r) 

(cs)n = -^^^iv - ^?o)diag(l, 0, g, 0), (14) 

where U is the energy density of string per unit length and fj. = qlJ is the tension. 

For current carrying strings and coarse-grained strings that are wiggling on microscopic scales, q generically takes a 
value different from unity. Here, we do not specify microscopic model of the string and will leave g as a free parameter. 

We begin by decomposing (cs)2^ii i^ito scalar, vector and tensor-type perturbations 

1 ^2A^ 1 



ics)T^J = yDiDj - ■^1^JD- j S(cs) + 27un(cs) + D.vj + DjVi + H,^, (15) 

where (cs)rij = gik(cs)T^j = a'^lik{cs)T), D^Vi = and H^j = 'f'^Hij = 0. Comparing Eqs. and (US]) then 

implies 

n(cs) = -^^e(7y-7?o), (16) 

^(CS) = ^ log r 6(77- 770), (17) 

V^ = 0, (18) 

Hrr^£^e{W0), (19) 

H,, = 0, (20) 

Hre = 0, (21) 

^^.. = -^e(,-,„), (22) 

H,e = 0, (23) 

Hee - ^^^e(r; - r;„) - ^e(r)e(r7 - 770). (24) 



Since (G'5)T^u is itself a first order quantity, {cs}Tj is simply given by 



(25) 



where H^j = ^'■^Hkj 



D. Energy-momentum tensor of a perfect fluid 



We assume that the perfect fluid obeys the equation of state given by P = c^p, where Cs is the sound speed. We 
then write the velocity perturbation as Ui = aDiV so that the perturbed four-velocity is given by 

(1 - $) 1 ■ \ 

^ '-,-D'v], u^ = (-a(l + $),aA«). (26) 

a a J 
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The nonzero components of the perturbed energy-momentum tensor are then 

STl = -dpr^, (27) 
={p + P)D,v, (28) 
5Tj=cl5p^5'^. (29) 

III. PERTURBATION EQUATIONS 

A. Scalar-type perturbations 

Six equations may be derived for the scalar-type variables, two from the conservation laws of the energy-momentum 
tensor and four from the Einstein equations. These constitute the basic equations that govern the evolution of the 
density contrast {5m = 5p/ p) and velocity perturbation {v) of the perfect fluid, the gravitational potential ($) and 
the curvature perturbation (vP). However, not all these equations are independent as they are related through Bianchi 
identities and it is possible to derive a closed differential equation for each variable. Before doing this, let us first take 
the following ansatz; 

<5m = <5,„e(7]- 770), w = we(77- ?7o), <!> = $6(77- 770), = *e(7/-77o). (30) 

This is reasonable since nothing must be excited before the phase-transition epoch, t] = rjo. The basic perturbation 
equations can be written in the following form; 

I + 3(1 + cl)ni> - cIdH = ^^'^^'^(^) , (31) 

$ + * = -2gGC/logr, (32) 

i?2$ = 2(l-g)GC/^ + ^A™, (33) 
r 2 

§ -'H^ = ^{l + cl)n^v, (34) 

where A,„ is defined by 

A,„ = ~S,n - 3(1 + cDnv. (35) 

Equation (|3T|) is a closed differential equation for ^. Once the time evolution of is determined, the other variables 
$, Am and v are detemined uniquely by Eqs. (f32|) - (f34| . The initial conditions for \P are given by 



^'(7]o,7') 0, '^{rio,r) = -2qGUHlogr. (36) 

and are obtained by requiring that the terms containing a delta- function S{ri — rio)i which appears in the perturbation 
equations after the change of variables Eg.ipn]). should vanish. 

We can simplify the master equation for ^ (PT|) by introducing a new variable, TZ, defined as 

7e = * + 2GC/logr. (37) 

The evolution equation for TZ is 

■R + 3{l + cl)nn-clD^TZ = 0,, (38) 

which is free from the source term given by the appearance of the string. Conveniently, all information regarding 
the string formation is contained solely in the initial conditions although, physically, this may have been expected 
since (at least within our current approximation), this produces a "source" term for the metric perturbations only 
instantaneously. The initial conditions for TZ are 

7^(77o,r) = 2GUlogr, TZ{qo,r) = -2qGUHlogr. (39) 
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B. Tensor-type perturbations 

As mentioned in Sec. II-B, the relevant evolution equation for the tensor perturbation is given by the r — r 
component of the Einstein equations, 

F + 2nF-F" --^F' ^^^Q{r,^r^o)- (40) 

It may be confirmed that the evolution equations for the other components can be derived from this. Taking the 
ansatz F(r],r) = F(r],r)Q{r] — 77), Eq. (j40p may then be rewritten as 

F + 2nF- F" --F' = ^^, (41) 
r 

and the initial conditions are, 

F(ryo,r) =i?(?7o,r) =0. (42) 



IV. SOLVING PERTURBATION EQUATIONS 

A. Scalar-type perturbations 

Let us try solving the master equation for (|38p subject to the initial conditions ([5^ . We first perform the 
Fourier-transformation for 72.; 

n^.r)= I -f^e^^'-^n.iv), (43) 

where x ~ r(cos 6, sin 6) is two-dimensional vector on the (r, 0)-plane. The evolution equation for each Fourier 
component is 

iik + '^^"^^^ Tlk + clk^-Rk - 0, , (44) 

■n 

where u is defined as v ~ j^^- In order to give the initial conditions for each 72,^, we must also Fourier-decompose 
the logr term in Eq. ()37p . Rewriting logr as 



\ogr= —^e'-^ak, (45) 

where ak are constants, we find that the identity log?- = yields Uk = — 27r/fc^. However, direct substitution of 
this into Eq. ()45p results in a divergence which does not depend on r. To extract r-dependcnt part wc must therefore 
modify the coefficients via the transformation ak ^ ock = — 27r/A;^~'^, where e is a positive infinitesimal. This choice 
of oik yields 

— 7^ = h logr -f- const. + 0(5). (46) 

Zi: k'^ ^ e 

so that the part which is independent of e gives the desired log r (the const, term is not important for our purpose). To 
find the desired solution for TZk, we keep e in the intermediate calculations before letting e -> in the final expression, 
thus picking up terms which are independent of e. With this prescription, the initial conditions for TZk can be written 

as 

AnGU ■ Anq^GU 
T^^ = -l^^ ^^ = ^237^- (47) 

It is then straightforward to derive the solution of Eq. (|44| , which is 



AttGUc if^^ 

Ti-kiv) = , " — {akju{_CskT]) -f bkn^{cskri)) , (48) 

rv Tj 
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where and 11^ are spherical Bessel functions of the first and second kind, respectively, and and 6^ are given by 

Ofe = -qiyn^{cskriQ) + Cskrion^+i{cskr]o), hu = qvju{csk-qn) - CskriQj^+i{cskr]o). (49) 

Using this solution with arbitrary values of the equation of state parameter Cs, and the string formation time 770, we 
can obtain the time evolution of TZiri, r) by integrating Eq. (j43|) over k numerically. 

However, let us consider a situation where an analytic evaluation of the integral (|43p may be obtained to good 
approximation. We begin by assuming that ?/o is sufficiently small, such that both r and the sound horizon distance, 
Cs?7, are much larger r]Q. This approximation applies to a comoving observer attached to the background FLRW 
universe (and whose position is therefore specified uniquely by the comoving radius r), at a time much later than the 
time of string formation. We do not assume a magnitude relation between CsT] and r, that is, CgT] can be larger or 
smaller than r. If Csi] > r the observer is inside the sound cone, whereas 0^77 < r means that the observer lies outside 
the sound horizon. 

Since the contribution to Ti.{r],r) from the part is suppressed by positive powers of 770, it vanishes in the limit 
770 — ^ 0. We therefore neglect the contribution from bk and keep only term. We then Taylor expand the spherical 
Bessel functions in Cskrjo and keep only the leading order term, so that 



{qv-2v-l)T{2v + l) 

2-r(. + i) ^^^'^^^ ■ 



With this approximation, 7?,(?y, r) becomes 



7e(?7, r) ~ 2^-''GU{qv ^2v-l 



T{2u + 1) 

'r(7. + i) 



dz 



(50) 



(51) 



where x is defined by x = r/{csr]). Using the formula for the integral of a product of two Bessel functions given in 
the appendix, we find that integration over z yields, 



dz z'^ 



^jQ{xz)jy{z) 



,.r(. + i)r(i-|) +2' 



for a; > 1, 



(52) 



and 



dz z 



e-u-l 



Jo{xz)ii,{z) 



V^2- 



^ - -1 
2 2) 



-F 



2' 



1 e 

2 2'' 



for X < 1, 



(53) 



where -F(a, &, c; x) is Gauss' hypergeometric function. After Taylor expanding these expressions and picking up terms 
independent of e by allowing e ^ (as before), we find that TZ is given by 



7^(7/, 7^ 



2{qv-2v- \)GU 



2v + l 

where we define the function Qu{x) as 



logr + const. 



e(a-- 1) + 



{qv-2v-l)GU 
2v + l 



(-21og77 + g,(x))e(l-x), (54) 



db 



-F b 



(55) 



6=0 



For later convenience, we here give an explicit expression for Q^, when = 1 in terms of elementary functions; 
Gi (x) = 2 fVl- x^ + " - + V^l - •^') - ^ + log 2V 



(56) 



This corresponds to the physically interesting case = 1/3, i. e. a radiation-dominated universe. Having determined 
TZ, we can immediately determine ^' from Eq. (|37p : 



2qvGU 
2v+l 



log 7' e{x -1) + GU 



21ogr + (^f_ij_il(-21og77 + g^{x)) 
2v -\-\ 



9(1 -x) 



(57) 



followed by $ from Eq. ((5^ . 
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B. Tensor-type perturbations 



By performing the transformation 'ip = rF ^ the evolution equation for F (j4ip may be rewritten in terms of the new 
variable "0 as 



y 2z/ ■ ,„ 1 1 , AqGU 



(58) 



A solution of this equation satisfying the initial conditions ipivoj''') = V'(^Oj'') = is 



■0(77,0 



dk Akir])Ji{kr), 



(59) 



where Ak (rj) is a function that obeys the following differential equation, 

Ak + —Ak + k^Ak = 4:qGU. 
Imposing the initial conditions Ak{rio) = Akirjo) = 0, we find that Ak{ri) is given by 



(60) 



Akir,) 



AqGU 



1 - krikrio^ ^ ) {3i,{krio)n^{kvi) - ji,.i{kri)ni,{kr]o)} 



Now, let us again consider the situation rj ^ riQ and take the limit 770 —> 0; Ak{ri) then becomes 



(61) 



AqGU r(2i. + l) j._i(fc77) 



(62) 



Correspondingly, F{r], r) can be written as 



AqGU f°° dz ( r(2i^+l) j^_i(z) 



^ 1 



y Jo z 



2T(t/- 



(63) 



where y is defined by y = r/ij. It is clear from the last equation that F depends on rj and r only through the ratio 
r/?7, which means that F exhibits self-similarity at late time Again, for later convenience, we here give F(y) for 
= 1 in terms of elementary functions; 



F{y) - AqGU 



^ + M^logi±^^ 



(65) 



To determine the evolution of the perturbation vari ables , Magueijo considered the closed differential equation for 
obtained from the continuity equation (V^T^ = 0) [l20j | 



A™ + (1 - ScDn^m - clD^Kn - -(1 + 2c2 - 3ct)n^A„, = 0. 



whose solution in the limit j] 3> ?7o and r ^ Cgrjo is given by 



A,„(77,r) 



AGU{~l + {q-2),^) ^ 



1-x 



(66) 



(67) 



^ One can arrive at F equivalent to Eq. l|63p by making a self-similar ansatz F{r], r) = F{y) from the outset and by solving the resultant 
ordinary differential equation for F given by 

3/2(1 _ y-^)F"(y) + y(3 + 2{u - l)y^)F'{y) + AqGU = 0. (64) 

Magueijo obtained the solution for the tensor-type perturbation by this approach [l2Cll |. 
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He set q = 1 throughout and ([67)) coincides with the solution given in |12Clj in this case. He then proceeded by solving 
Eqs. ([5^ and p3p for the physically interesting case oi i/ = 1 (i.e. for a radiation-dominated Universe), to obtain 
expressions for $ and 4*. However, the expression obtained in |l20l | does not satisfy the master evolution equation for 
^, Eq. ([nj), obtained from the trace part of the Einstein equations. The inconsistency can be resolved by adding a 
time dependent function which is a homogeneous solution of the Poisson Eqn. (|33p . The general solution of Eq. ([33)) 
which is continuous and whose first derivative with respect to r is continuous at r = CsTj is 



GU 



4 2 
Qiv) - 3 log 2:+ - 



2\ogx + - - 21og2 



9(1 - a;) 



(68) 



where Q{r]) is an arbitrary function. For the specific choice Q{r]) = '■, 
in jl20l |. but the corresponding expression for ^, calculated from Eq 
of the following differential equation; 



77 



(3 log 2 — 4), this solution reduces to that given 
([32)) . satisfies Eq. ([3T)) only if Q{ri) is a solution 

(69) 



The general solution of Eq. 



691) is given by 

Qiv) 



-g (log'? + 3 - 10g2 



(70) 



and substituting this back into Eq. ()68p yields an expression for <!> which coincides with the one obtained previously 
by directly solving Eq. ([5T)) . 



V. PROPERTIES OF THE PERTURBED SPACETIME 

The solutions obtained in the previous section now allow us to study the properties of the metric perturbations 
excited by the string formation. For simplicity, throughout this section, we confine our analysis to the case of a 
radiation-dominated universe, i. e. u ~ 1, in which case analytic calculations are feasible to some extent. However, 
we expect that the qualitative picture we obtain will hold generally true in the presence of background fluids governed 
by an equation of state, P = c^p, with 7^ 1/3 and, perhaps, even in more exotic cases. 

The tensor-type and curvature-type metric perturbations given by Eqs. ()57p and ()65p contain step functions which 
are discontinuous at the sound horizon r = CgT] and the cosmological horizon r = rj, which implies that the properties 
of the perturbed spacetime differ qualitatively within each region. We therefore split the spacetime into three regimes, 
I) outside the cosmological horizon r > rj, 11) between the cosmological and the sound horizons c^ry < r < j] and HI) 
inside the sound horizon r < Csf], and consider the properties of the perturbed spacetime for each regime separately. 



A. (I) Outside the cosmological horizon 

In this regime, the perturbed metric is given by 

(1 + 2$/) drf + {I + <l>i + Fi) dr^ + (l + - Fi) dO^ (1 -(- $/) dz^ 



ds^ = rf 

where $/ and Fj are defined according to 



AqGU , ^ 2qGUr]^ 

■ log r, Fi = 



3 3r2 • 

However, performing the coordinate transformation (77, r) — > {ri{T,R), r(T, i?)), where 

V = T + —^TlogR, r^R+ , 

we find that this metric reduces to the unperturbed FLRW metric, 

ds^ = {-dT^ + dR^ + R^dO^ + dz^) . 



(71) 

(72) 

(73) 
(74) 



* Since the solutions are correct only up to first order in GU , we must keep in mind that we neglect any second or higher order terms in 
GU in the new metric. 
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It can be confirmed that the density perturbation and velocity perturbation of the perfect fluid also vanish in the new 
coordinates. Therefore, despite the initially highly noncausal appearance of the perturbation variables ([7^ . we see 
that information regarding the string formation does not travel beyond the light cone, and the FLRW metric remains 
valid beyond the causal horizon. One may wonder how it is possible for the gauge-invariant scalar and tensor-type 
perturbations, which were nonvanishing in the original coordinates, to give vanishing perturbations beyond the light 
cone in the new coordinates. In fact, it is true that, even after the coordinate transformation, neither the scalar-type 
perturbations, nor the tensor-type perturbations vanish individually. However, in the new coordinates they are exactly 
equa l in magnitude, but opposite in sign and thus cancel, yielding vanishing net perturbations. As pointed out in 
[l2C| , the apparent violation of causality originates from a fact that the geometrical splitting of a perturbation into 
scalar and tensor type is a nonlocal operation and manifest causality is recovered after all types of perturbations are 
taken into account. 



B. (II) Between the sound and the cosmological horizon 



In the middle regime, the metric can be written as 



I 2 2 

as = 1] 



(1 + 2$/) dif + (1 + $j + Fii) dr^ +r^{l + ^i + Fn + rF'^) d9^ + (1 + - 2Fii - rF'j^) dz"' 



where ' indicates the derivative with respect to r and Fn is defined by 

1 



F, 



II 



iqGU 



1 1/3, 



r 

y = -■ 



, (75) 



(76) 



Since the matter excitation should be absent outside the sound horizon, we expect that only gravitational degrees of 
freedom, i. e. tensor-type perturbations, exist in this regime. Given that all information regarding the tensor-type 
perturbations is contained within the spatial components of the metric, it would be convenient to perform a coordinate 
transformation for which the resultant, new metric, has only spatial perturbations. In fact, this is already achieved 
by the coordinate transformation given in Eq. (j73p , which yields 



ds^ 



dT^ + (1 + gii)dR^ + R^{1 + Gii + Rg'ii)de^ + (1 - 2^/7 - Rg'ij)dz' 



where Qn is defined as 



gn = -GUq 



2^2 



1) + log K — 3/t2 log I 



3k2 



(77) 



(78) 



and K, = ^ for the metric perturbations. The density perturbation and velocity perturbation of the perfect fluid 
still vanish in the new coordinates. As expected, the metric perturbations in the new coordinates still satisfy the 
transverse and traceless conditions, showing explicitly that only tensor degrees of freedom are excited in this regime. 
It is interesting, however, that they depend on T and R only through the ratio k, which means that the gravitational 
waves caused by string for mati on exhibit self-similarity, i.e. waves at (T, R) look the same as those at (6T, hR), where 
& is a numerical constant Il26| . Magueijo assumed a self-similar ansatz in this regime and our results confirm that 
this is correct. However, [12Cll | also states that the gravitational potential in the region c^ry < r < is that of a 
static repulsive rod, whereas the metric ([77)) shows that a potential is not present between the sound horizon and 
the cosmological horizon. Now let us consider the embedding of the two-dimensional surface spanned by the R and 6 
coordinates at fixed T and z. The line element of the subspace is given by 



Introducing a new radial coordinate i?. 



wc find that this can be rewritten as 



(1 + gH)dR^ + R^{i + gii + Rg'n)de'^ 



R^R[l + \gn + \g'u ) , 



(1 + 4qGC/Vl - ~^'^)dR? + R^dQ'^ 



. R 

K = —■ 

T 



(79) 



(80) 



(81) 
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R/T 



Figure 1: Plot of h{K) outside the sound horizon using the normahzation yJqGUT = 1. 



Therefore, 2ttR is the circumference of a comoving circle whose radial coordinate is specified by R. The height of 
the embedding surface in the three-dimensional space, relative to that for the region outside the light cone, which we 
denote h{R), is then given by 



h{R) = 2^qGUT 



V^r(|) R^ / 1 1 3.^^ 

2r(|) T \ 4'2'2'T2 



(82) 



We find that the first derivative of h{R) aX T = R vanishes but the second derivative diverges. A plot of h{R) is 
given in Fig. [T] Next, let us consider how the infinitesimal distance between two particles following initially parallel 
geodesies changes as they across the gravitational wave front at the c osmo logical horizon. The change of distance is 
described by the geodesic deviation equation (sec, for example, Wald [l27l |). 



dr2 



(83) 



where is the infinitesimal displacement vector between the two particles and r is the proper time.^ Assuming that 
the particles are initially attached to the comoving coordinates when they are outside the light cone, the component 
R' 

ooi gives the dominant force for the displacement along i-th axis. The relevant components are therefore. 



00 



009 



ooz 



2qGU (k^ + 2) vT 



2qGU (2k4 



1 

1 

1 2qGUK 

r2 



(84) 
(85) 
(86) 



The first term, — l/T^, common to all the components above is merely the tidal force due to the cosmological expansion. 
It is simply the second time derivative of the scale factor and the minus sign is due to decelerating expansion of the 
Universe. The second term in each component represents the additional tidal force induced by the gravitational waves. 
We find that, at the cosmological horizon where k = 1, the i?-component of the Riemann tensor vanishes while the 
9 and Z-components diverge toward positive and negative infinity, respectively. There is therefore a gravitational 



5 The Riemann tensor is defined by R"^^^ = d^Ff^^ - d.^F"^^ + T^^^F^^^ - T^^^^F^^^. 
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Figure 2: The motion of free particles after crossing the cosmological horizon. 



shock front at the edge of the hght cone and, although the tidal force diverges, this does not imply that two particles 
separated by a small distance undergo a large displacement as they cross the cosmological horizon. The divergence 
is weak enough to make the solution of the geodesic deviation equation, i.e. the double time-integration of the tidal 
force, finite. The motion of a ring of free particles hit by the gravitational shock wave at the cosmological horizon is 
schematically plotted in Fig. [5] The reason that the ring is elongated along the angular direction as the particles cross 
the cosmological horizon is intuitively understood as follows. Let us consider a ring which consists of free particles 
with the same radial coordinate Rq, separated by equal angles, and vary Rq from Rq > T to Rq < T. This procedure 
represents how the ring is distorted as it crosses the cosmological horizon. Before the ring hits the cosmological 
horizon, the ratio of the rate of decease of the circumference of the ring to that of the radius is given by 27r and the 
particles move closer together at a constant rate. After the ring crosses the horizon, due to the distortion of space 
caused by the gravitational wave, the rate of decrease of the physical radius (which is the length of the solid line in 
Fig- nil becomes greater than that of the circumference. This can be understood by considering the limiting case 
where the slope of the curve for R/T in Fig. [T] is very steep. In this case, even a significant translation along the 
solid line results in a tiny shift in R. Therefore, the rate at which particles on the ring move closer to each other is 
slower inside the horizon, than it is outside, which is equivalent to saying that the particles have been affected by the 
outward tidal force that stretches the distance between two particles. 

C. (Ill) Inside the sound horizon 

In the final regime we consider the region of spacctime inside the sound cone, in which the metric can be written as 



ds^ = if 



- (1 + 2$///) drf + (1 + 2^-///) { (1 + Fii) dr^ + r^l + Fn + rF'^) d0^ + (1 - 2F// - rF^j) dz^} 



, (87) 



where ^/// is defined as 



= GU 



21ogr + ^^^(-21og7y + 5i(x)) 



and is then given by Eq. ((32|). We find that the following coordinate transformation, 

3GC/ - - 1)Z2) GU , , , 
V = T ^ ^ ^__(^_3)r(21ogr-l), (89) 

GU 

r = R + GU(q + 2)RlogR i?(2g(2 + 31og2) + ISlogT + 9), (90) 

6 

GU 

z = Z+— Z(3-13g + 6glog4 + 18(g-l)logr), (91) 
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yields a new metric which reduces to 



{ - l + 'iGU{q-\)\ogR]dT^ + dR^ + R^{l-AGU{\ + q)]de^ + {l + AGU{q- I) log R]dZ^ 



(92) 



in the subhorizon hmit, R/T <^ 1. This is the metric for an eternal cosmic string and coincides with the one given by 
Vilenkin, Gott and Hiscock in previous studies (53l - [55| . Therefore, as expected, the spacetime sufficiently close to a 
noneternal cosmic string resembles that for the eternal string. The nonvanishing time-time component of the metric 
perturbation leads to a gravitational acceleration for test particles given by _?:£lIJ^ii3l ^ \i jg therefore an attractive 
force for g < 1 and becomes repulsive for q > 1, which also makes sense intuitively since < g <C 1 corresponds to 
a nonrelativistic string and <? > 1 corresponds to a tension-dominated string. For the Nambu-Goto string, for which 
(7=1, there is no gravitational force. In all cases, there is a deficit angle around the string of 47rGC/(l -I- q). 

The general form of the metric in the new coordinate system, without the subhorizon approximation, is very 
complicated. For simplicity, we here consider only the Nambu-Goto string, for which q = 1. In this case, the metric 
perturbations, which are valid in the whole region inside the sound horizon, can be written as 



T-TT, (777) 



h 



RR,{III) 



■XIII) 



^ZZ,(III) 



9 

GU 
9^ 



GU ( Vl - x^x^ - 4^1 - 



31og( v/l-a:2 + lj + 4- log 

2-41 + 271og3 + 421og2 



(93) 



6^9 - 3x2 _ (^4^9 _ 32;2 + 32 
+Zx^ {i log (-x^ + 2v/9 - 3x2 + 6^+8 log (^-Jl-x"^ + l) + 6 log + (8^1-2:2 -Q^x^ + 18 
6 f V9 - 3.x2 - 3) -x^ f 8v/9 - 3a;2 -f 32 ^1 - a;2 + 13 + 9 log 3 + 42 log 2 
2 + 6) +81ogfv/l-a;2 + lll + f8v/l-x2-9 



GU 
9^ 



3x^ (31og (-x^+2\/9 - 3x 



(94) 



(95) 



-GU 
9 



(sVl-x^ -9^x^ + 2 (eV^ - 3a;2 - IGv^l - a:2 - 2 + 6 log 2 + 9 log 3 
18 log (^-x^ + 2^9 - 3a;2 6^ 24 log (^\/l-x^ + 1 



(96) 



and all other components vanish. It is interesting to see that all the components of the metric perturbation depend on 
T and R only through x. Therefore, as in the case for the middle regime, the perturbations also exhibit self-similarity 
inside the sound cone. Unlike when we approach the edge of the light cone, we find that the Riemann tensor does not 
diverge at the sound horizon and the relevant components are given by 



^ OOR 



1 



2{15 + 7VQ)GU 



R', 



rp2 Qrp2 ' OOe 

Sufficiently close to the string core, these reduce to 



1 , (17V6-30)G[/ z 
9T2 ' 



R 



OOR 



1 



6GU 

J<2 



+ 0{x'), R 



R 

OOR. 



1 

T2 



QGU 



+ 0{x^), R 



z 

ooz 



)_ _ {IQ + V&)GU 
r2 3T2 • 



= -- + 0{x^). 



(97) 



(98) 



There therefore remains a residual outward tidal force along both the radial and angular directions, which does not 
vanish even in the limit a; — >■ but which, of course, still vanishes in the eternal limit, T — >■ 00. 



VI. CONCLUSION 



We derived the analytic form of the metric perturbations on an FLRW background, excited by the formation of a 
straight cosmic string under the wire approximation. Our results are valid for a general string with differing tension 
and energy per unit length. For the Nambu-Goto string, our metric does not coincide with that given in previous 
work due to a differing choic e of boundary conditions for the Poisson equation. We found that, as a result of this 
choice, the solutions given in |l2C| for the gravitational potential and curvature perturbation variables do not satisfy 
the trace part of the Einstein equations. 

By performing the appropriate coordinate transformation, we explicitly verified that the spacetime outside the light 
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cone remains unperturbed and governed by the FLRW metric. Whereas it is nontrivial to demonstrate nonviolation of 
causality in the original metric, it is manifestly respected in the new coordinate system. At the cosmological horizon, 
we found that there exists a gravitational shock giving rise to a divergent tidal force along both the angular direction 
and the direction of the string axis, z. The shock distorts a small ring of test particles placed in the 9 — z plane by 
stretching it along 0-direction and shrinking it along the z-direction. The degree of divergence is weak, so that the 
resultant distortion remains finite. Inside the cosmological horizon, the metric perturbations exhibit self-similarity at 
late time, i.e. they depend only on a ratio of the radial coordinate to the conformal time. This is reasonable since 
the Hubble distance is the only characteristic length scale in the system. Deep inside the horizon, we found that the 
metric tends to that for an eternal string at late times, which has been used to search for the signatures of cosmic 
strings in astronomical observations. 
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Appendix A: Formula for integral of a product of two Bessel functions 

In this appendix, we provide a formula for the integral / defined by 

/•OO 

/ = / x^J^{ax)J^{bx). (Al) 

For a > > 0, we have 

T — - 

a>^+^+^T{v+\)T{\ -13) V a 
where F is the usual Gauss hypergeometric function and a and /3 are defined by 

a^^{X + H + u + /3= i(A-^ + !/ + l). (A3) 







For a = 6 > 0, we have 



and for 6 > a > 0, 



I . 2^rian-X) 

a^+^r{l ~ p)r{iy - a + l)riiy - p + 1) ^ ' 



-F\ci,a-v,ii + \--r\. (A5) 



5A+M+ir(^+ l)r(:/-a + 1) V' "^ 62 
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